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ABSTRACT 

The process of baryon number transfer due to string junction propagation in rapidity 
is considered. It has a significant effect in the net baryon production in pp collisions 
at mid-rapidities and an even larger effect in the forward hemisphere in the cases of 
np and 7p interactions. The results of numerical calculations in the framework of the 
Quark-Gluon String model are in reasonable agreement with the data. 
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1. INTRODUCTION 



The Quark-Gluon String Model (QGSM) and the Dual Parton Model (DPM) are 
based on the Dual Topological Unitarization (DTU) and describe quite reasonably many 
features of high energy production processes, including the inclusive spectra of different 
secondary hadrons, their multiplicities, KNO-distributions, etc., both in hadron-nucleon 
and hadron-nucleus collisions [l], 0, [|. High energy interactions are considered as pro- 
ceeding via the exchange of one or several pomerons and all elastic and inelastic processes 
result from cutting through or between pomerons ||. The possibility of exchanging a 
different number of pomerons introduces absorptive corrections to the cross sections 
which are in agreement with the experimental data on production of hadrons consisting 
of light quarks. Inclusive spectra of hadrons are related to the corresponding fragmenta- 
tion functions of quarks and diquarks, which are constructed using the reggeon counting 
rules ||. 

In the present paper we discuss the processes connected with the transfer of baryon 
charge over long rapidity distances. In the string models baryons are considered as 
configurations consisting of three strings attached to three valence quarks and connected 
in a point called "string junction" |7|, ||. Such a configuration corresponds in QCD to 
the following gauge-invariant operator 

^(xi)^(x 2 )^ fc (x 3 )G[P(x 1 , x%G[P(x 2 , x)]i,G[P(x 3 , , (1) 

where 



G[P{x u x% 



Texp{g f A^x)dx^ 

JP(xi,x) 



(2) 

J J 

It is very important to understand the role of the string-junction in the dynamics of 
high-energy hadronic interactions. Now we have several different experimental results 
concerning such processes. First of all the data clearly show that in the forward 
hemisphere the number of secondary protons produced in n + p interactions is significantly 
larger than the number of p produced in ir~p collisions. This difference can not be 
described M without the assumption that the baryon charge is transferred from the 
target proton to the pion hemisphere. 

Similar data on the differences of p — p yields in |(vr + j9 + n~p) collisions at 158 GeV/c 
were presented by the NA49 Coll. |T0|. 

A second group of data concerns the energy dependence of the differences in yields of 
the protons and antiprotons at 90° (i.e. at xf = 0) at ISR energies [fPT| . 

Another sample of data includes the measurements of hyperon production asymme- 



tries in 500 GeV/c tt -nucleus interactions |12 



Finally, the proton-antiproton asymmetry in photoproduction was recently measured 
at HERA 0. 



In this paper we present a simultaneous description of all these data and extract 
information on the properties of the string junction dynamics. 
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2. INCLUSIVE SPECTRA OF SECONDARY HADRONS IN QGSM 



As mentioned above high energy hadron-nucleon and hadron-nucleus interactions are 
considered in the QGSM and in DPM as proceeding via the exchange of one or several 
pomerons. Each pomeron corresponds to a cylindrical diagram (see Fig. la), and thus, 
when cutting a pomeron two showers of secondaries are produced (Fig. lb). The inclusive 
spectrum of secondaries is determined by the convolution of diquark, valence and sea 
quark distributions u(x, n) in the incident particles and the fragmentation functions G(z) 
of quarks and diquarks into secondary hadrons. The diquark and quark distribution 
functions depend on the number n of cut pomerons in the considered diagram. In the 
following we use the formalism of QGSM. In the case of a nucleon target the inclusive 
spectrum of a secondary hadron h has the form M: 



where the functions 0^(x) determine the contribution of diagrams with n cut pomerons 
and w n is the probability of this process. Here we neglect the contributions of diffraction 
dissociation processes which are comparatively small in most of the processes considered 
below. It can be accounted for separately |I|, 0, . 
For pp collisions 



<Wc) =fi q (x + ,n)f^x.,n)+f^x + ,n)f^x.,n)+2(n-l)f^x + ,n)f^x^n) , (4) 



where f qq , f q and f s correspond to the contributions of diquarks, valence and sea quarks 
respectively. They are determined by the convolution of the diquark and quark distri- 
butions with the fragmentation functions, e.g., 



In the case of a meson beam the diquark contributions in Eq. (4) should be changed by 
the contribution of valence antiquarks: 



= f^x + ,n)f^x.,n) + f h q {x^n)f h qq {x^n) + 2(n - l)fi{x+, «)/*(*_, n) . (7) 



The diquark and quark distributions as well as the fragmentation functions are deter- 
mined from Regge intercepts. Their expressions are given in Appendix 1. 

The net baryon charge can be obtained from the fragmentation of the diquark giving 
rise to a leading baryon (Fig. 2a). A second possibility is to produce a (leading) meson 
in the first break-up of the string and the baryon in a subsequent break-up (Fig. 2b). 




(3) 




(5) 
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As discussed above, in the approach ]7|, ^ the baryon consists of three valence quarks 
together with string junction (SJ), which is conserved during the interaction^]. 

This gives a third possibility for secondary net baryon production in non-diffractive 
hadron-nucleon interactions (Fig. 2c). 

The secondary baryon consists of the SJ together with two valence and one sea quarks 
(Fig. 2a), one valence and two sea quarks (Fig. 2b) or three sea quarks (Fig. 2c). The 
fraction of the incident baryon energy carried by the secondary baryon decreases from 
a) to c), whereas the mean rapidity gap between the incident and secondary baryon 
increases. 

The probability to find a comparatively slow SJ in the case of Fig. 2c can be estimated 
from the data on pp annihilation into mesons (see Figs, lc, d). This probability is known 
experimentally only at comparatively small energies where it is proportional to s asj_1 
with asj ~ 0.5. 



However, it has been argued [OJ] that the annihilation cross section contains a small 
piece which is independent of s and thus asj ~ 1. Irrespectively of the value of asj, the 
contribution of the graph in Fig. 2c corresponds to annihilation and, thus, has a small 
coefficient which will be denoted by e. In our calculation we shall use 

a S j = 0.5 (8) 

and will treat e as a free parameter. A contribution with asj ~ 1, if present at all, is 
expected to be important only at very high energies (see a discussion of this point in 
section 5). The values of asj and e can only be determined with the help of accurate 
data. 



3. COMPARISON WITH THE DATA 

The mechanism of the baryon charge transfer via SJ without valence quarks (Fig. 2c) 
was not accounted for in previous papers U pj, f||. 

The data at comparatively low energies (y/s ~ 15 4- 40 GeV) can be described with 
a sj — 0.5. Unfortunately, the value of e cannot be determined in a unique way. Many 
sets of data can be described with e = 0.05. However, other data sets favor a value four 
times larger, e = 0.2. Because of that, we will present our results for these two values 
of e, leaving all other parameters in the model unchanged. Thus, the results for meson 
and antibaryon production are the same in the two cases. 

The inclusive spectra of secondary protons and antiprotons produced in pp collisions 
at lab. energies 100 and 175 GeV || are shown in Figs. 3a and 3b together with the 
curves calculated in the QGSM. The agreement of our results with e = 0.05 with the data 
is quite reasonable, on the same level (or even better) as in the previous papers which 

'At very high energies one or even several SJ pairs can be produced. 
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did not incorporate the SJ mechanism of Fig. 2c. The variant with e = 0.2 gives too 
large multiplicity of secondary protons at small Xp. The description of the antiproton 
yields is reasonable. 

The data of Ref. [15] are in some regions of xp in disagreement with the data of Ref. 
H as well as with our calculations. However, one can see from Figs. 3c and 3d that 
the experimental points can not be described by any smooth curve. Probably this is 
connected with the use of different detectors for different xp regions and to the influence 
of the trigger in |I5] . Again, the value e = 0.05 is preferable. 

The data on secondary proton and antiproton production in pp collisions at ISR 



energies [FT] at 90° in c.m.s. are presented in Figs. 3e and 3f. Their differences, which 
are more sensitive to the baryon charge transfer, are presented in Fig. 3g. One can 
see that the last data, as well as the yields of protons and antiprotons separately, are 
described quite reasonably by QGSM with e = 0.2. However, it is necessary to note, 
that the systematic errors in [O] are of the order of 30 %, so the value e = 0.05 can not 
be excluded. Thus the disagreement between ISR data [|TTJ and more recent data || [15 



on spectra of protons in pp-collisions does not allow to determine uniquely the value of 

e. 

The data on baryon production in the pion fragmentation region || are presented 
in Fig. 4. The spectra of antiprotons produced in i\~p collisions, shown in Fig. 4a, 
allows one to fix the fragmentation function of a quark into baryon/antibaryon. If the 
contribution of the baryon charge transfer were negligibly small, the inclusive spectra 
of reactions ir^p — > pX and 7r + p — > pX in the pion fragmentation region would be 
practically the same M. Actually, the data for the second reaction are significantly 
higher than for the first one providing evidence for the baryon charge transfer due to the 
SJ propagation. The difference of the inclusive spectra in the two considered processes 
allows one to estimate quantitatively the contribution of the baryon charge transfer, 
and the parametrization of these processes given in Appendix leads to a reasonable 
description of proton yields in ir + p collisions with e = 0.2. The value e = 0.05 seems 
to be too small here. Note, however, that our factorized formulae, Eqs. 4 and 7, imply 
that the slowed down proton is made out of SJ and three sea quarks (see Fig. 2c). Thus, 
they should be modified for the reaction ir + p — > pX in the pion fragmentation region due 
to the possibility of SJ recombination with a pion valence quark - which would change 
the ratio of protons and neutrons. A simple estimate based on the quark combinatorics 
given in Appendix 1 leads to an increase of the proton yield in the reaction n + p — > pX 
by about 50 %. In this case the value of e = 0.05 can be consistent with the data. 
This problem disappears if we consider the sum of n + p and 7t~p. Preliminary data on 
\{k + P + vr~p) — > (p — p)X were obtained by the NA49 Coll. | I0| . One can see from 
Fig.4c, that they are in good agreement with our calculations with e = 0.05 and in total 
disagreement with e = 0.2. 

The data on A and A production in pp collisions |U| [17], 18[ [19J, presented in Fig. 5, 
are also in agreement with QGSM and the value e = 0.05 is again favored. 



5 



In Fig. 6 we show the data [EJ on the asymmetry of strange baryons produced in ir 
interactions^ at 500 GeV/c. The asymmetry is determined as 

for each xp bin. 

The theoretical curves for the data on all asymmetries calculated with e = 0.05 are 
in reasonable agreement with the data with the exception of at positive xp. 

However in this region the effect of recombination of SJ with valence (d) quark of 7r~ 
can be important and can lead to an increase of cross section analogous to the one 
in 7i + p — > pX discussed above. In the case of Q/fl production we predict a non-zero 
asymmetry in agreement with experimental data. Let us note that the last asymmetry is 
absent, say, in the naive quark model because Q and Cl have no common valence quarks 
with the incident particles. The value e = 0.2 seems to be excluded. 

Preliminary data on p/p asymmetry in ep collisions at HERA were presented by the 



HI Collaboration [13]. Here the asymmetry is defined as 



A ° = 2 t^t • (10) 

i.e. with an additional factor 2 in comparison with Eq. (9). The experimental value of 
Ab is equal to 8.0 ± 1.0 ± 2.5 % [|1| for secondary baryons produced at xp ~ 0.04 in the 
7P c.m. frame. QGSM with e = 0.05 predicts here 2.9 %, i.e. a smaller value, whereas 
the calculation with e = 0.2 gives the value 7.7 %, in agreement with the data. 



5. CONCLUSIONS 

We have demonstrated that experimental data on high-energy hadronic interactions 
support the possibility of baryon charge transfer over large rapidity distances. Probably, 



the most important are the baryon-antibaryon asymmetry in and [[L2] in n~p colli- 
sions, where both secondary particles have no common valence quarks with the incident 
particles. This asymmetry is provided by baryon charge transfer due to string junction 
diffusion. 

Also the production of net baryons in up and 7p interactions in the projectile hemi- 
sphere provides good evidence for such a mechanism. 

As for the values of the parameters a sj and e which govern the baryon charge transfer 
we have seen that the data, at comparative low energies (y/s ~ 15 -i- 40 GeV), without 
providing a clear cut answer, strongly favor the values asj = 0.5 and e = 0.05. Indeed, 



"'"These data were obtained from pion interactions on a nuclear target where different materials were 
used in a very complicated geometry. We assume that the nuclear effects are small in the asymmetry 
ratio (9), and compare the pion-nucleus data with calculations for n~p collisions. 
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as discussed above, the data which favor e = 0.2 have either a large systematic error |TT 
or correspond to situations where effects, not incorporated in the present version of the 
model, are expected to be important (as in n + p — > pX ; see the discussion in Section 3). 

At HERA energies the situation is different. Here the observed asymmetry is better 
described with e = 0.2. Note, however, that HERA data are preliminary and have 
rather large errors (the prediction for asymmetry with e = 0.05 deviates from the data 
by about two standard deviations). If this discrepancy is real can it be solved ? The data 
a t y/s ~ 40 GeV that we have examined provide information for the transfer of baryon 
charge over about five rapidity units. In contrast, the HERA data provide information 
for the corresponding transfer over more than seven units. In ref. [^TJ HERA data were 
described under the assumption that there is a component with ctsj ~ 1 and a very 
small strength. Another possibility would be to use asj — 0.5 and to assume, following 



Ref. [p2|-[24 that the strength of the SJ transfer contribution increases with the number 
of inelastic collisions. In this way the asymmetry would increase with energy since the 
average number of inelastic collisions, determined from Eqs. (42) and (43), increases 
with s. (This mechanism also allows to explain the increase of stopping in heavy ion 



collisions with increasing centrality) p2|-[26|. 



In conclusion, the study of the transfer of baryon charge over large rapidity distances 
is very important. Good experimental data in pp, pA and AB collisions for different 
centralities are needed in order to understand the dynamics of this mechanism. 

Acknowledgments: We thank H.G.Fisher from NA49 Collaboratiomn for providing 
the preliminary data in Fig. 4c prior to publication. 

This paper was supported by grant NATO OUTR. LG 971390 and NATO PSTCLG 
977275. 



APPENDIX. QUARK AND DIQUARK DISTRIBUTIONS AND THEIR 
FRAGMENTATION FUNCTIONS 

In the present calculations we use quark and diquark distributions in the proton for 
one-pomeron exchange given by the correspondent Regge behaviour [|]] 

u uu{ x j' n ) = C uu X aR 2a - B + 1 (l — x) aR , (11) 

u ud {x,n) = C ud x aR - 2aB (l-x)- a * , (12) 
u u (x,n) = C u x~ aR (l- x) aR ~ 2aB+1 , (13) 
u d (x,n) = C d x- aR (l-x) n+aR - 2aB , (14) 

(15) 

In the case of multipomeron exchange the distributions of valence quarks and diquarks 
are softened due to the appearance of sea quark contribution. Aslo it is necessary to 
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account that d-quark distribution is more soft in comparison with w-quark one. There 
is some freedom || how to account for these effects. We use the simplest way and write 
diquark and quark distributions (for cur = 0.5 and «b = —0.5) as 

u uu (x,n) = C uu x aR - 2aB+1 {l-x) fracA3 ( n - 1) - aR , (16) 
u ud (x,n) = C ud x aR - 2aB (l-x) n ~ l - a * , (17) 

(18) 

u u (x,n) = C u x- aR (l-x) n+aR - 2aB+1 , (19) 

(20) 

= C d x- aR {l-x)l {n - 1)+aR - 2aB+1 , (21) 

(22) 



u d \x,n 



Uu(x, n) = Uj(x, n) = C^x aR 

x [(1 - x )"+««-2°b-i _ tf / 2 (l - x y+** R -2a B -i] , n > i ; (23) 

u s (x,n) = C s x-° lR {l-x) n+2aR ~ 2aB ~ 1 , n> 1 . (24) 
In the case of a pion beam we use 

u q {x,n) =C q x~ aR {l-x) n - aR - 1 ,Ug{x,n) =C gX - aR {l-x) n - aR - 1 , (25) 
for valence quarks, and 



u(x, n) = Cx~ aR (l - x) n - aR - l [l - 5VT^~x\ ,n>\, 

u s {x,n) = C s x~ aR (l -x) 71 - 1 , n > 1 , (26) 

for sea quarks, where 5 = 0.2 is the relative probability to find a strange quark in the 
sea. The values of aR and as are given in Table 1. The factors C{ are determined from 
the normalization condition 

/ uAx, n)dx = 1 . (27) 
Jo 

and sum rule 



/ 52 u i( x -> n)xdx = 1 . (2£ 

JO ; 



is fulfilled. 

The fragmentation functions of quarks and diquarks were changed a little in compari- 
son with Refs. ]T|, |4[] to obtain a better agreement with the existing experimental data, 
because the old functions JTJ, |j] correspond to £ = 0. We use the quark fragmentation 
functions in the form 

Gl = G p d = Ml - z)*+°*-**»(l + a lZ 2 ) , Gl = G p d = (1 - ^)G£ (29) 
Gt,= G d = —(1 - z) A "G£ , Gj = (1 - s)G$ (30) 

"TV 
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Gf = ^(1 - z) Aa Gl , Gf = Gl = (1 - *)Gf (31) 

G n a = G? = G° = G? = ^(1 - z)^Gl , (32) 

with 

Aa = a p - = 1/2 , A = 2a' < p 2 >= 0.5 . (33) 

Diquark fragmentation functions have more complicate forms. They contain two 
contributions. The first one corresponds to the central production of a BB pair and can 
be described by the previous formulas. They have the form: 

G uu = G ld = GL = G p ud = Ml - z )A--*h<i-«b) , (34) 



z^A _ /^A _ /oA _ ^A _ A /-i _ "\ Aa/ip /or\ 

"AT 

G„ M = = = Gud — — (1 ~~ ^) Aa G«« 5 (36) 

with the same Aa Eq. (29). 

The second contribution is connected with the direct fragmentation of the initial 
baryon into the secondary one with conservation of the string junction. As discussed 
above, there exists three different types of such contributions (Figs. 2a-2c). Obviously, 
in the case of H production only two possibilities exist with string junction plus either 
one valence quark and two sea quarks or three sea quarks. In the case of production of 
a secondary baryon having no common quarks with the incident nucleons only the bare 
string junction without valence quarks can contribute (Fig. 2c). 

All these contributions are determined by Eqs. similar to Eq. (6) with the corre- 
sponding fragmentation functions given by 

G p uu = G p ud = a NV G[v e(l - z) 2 + v q z 3 ' 2 {l - z) + v qq z 2 } , (38) 

Gtd = aW^Ml - z ? + ^ 3/2 (! -*) + V'K 1 " *) Aa , G t = (1 - z)Gt d (39) 

Gfsj = a N V~z[v Q e{l - zf + v q z 3 ' 2 {l - z)]{\ - z) 2Aa , G% SJ = (1 - z)G%gj , (40) 

G% = a N v e^(l - z) 2+3Aa . (41) 

The factor y/z is really ^ 1_asj with asj = 1/2 (8). As for the factor yfzz 3 / 2 of the 
second term it is just 2(ctR — as) [0- For the third term we have just added an extra 
factor z x l 2 . 
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The probabilities of transition into the secondary baryon of S J without valence quarks, 
I 3 , SJ plus one valence quark, J 2 , and SJ plus a valence di quark, 7 l5 were taken from 
the simplest quark combinatorics [21|. Assuming that the strange quark suppression is 
the same in all these cases, we obtain for the relative yields of different baryons from SJ 
fragmentation without valence quarks : 

I 3 = 4L 3 : 4L 3 : 12L 2 S : 3LS 2 : 3LS 2 : S 3 (42) 

for secondary p, n, A + E, H°, S~ and Q, respectively. 
For 1 2 we obtain 



I 2u = 3L 2 : L 2 : ALS : S 2 : (43) 



and 



7 M = L 2 : 3L 2 : 4LS : : S 2 (44) 

for secondary p, n, A + E, S° and 
For Ji we have 

and 



J luu = 2L : : S (45) 



J lud = L : L : S (46) 

for secondary p, n and A + E. The ratio S/L determines the strange suppression factor 
and 2L + S = 1. In the numerical calculations we used S/L = 0.2. 

In agreement with experimental data we assume that E + + E~ = 0.6A |23j] in Eqs. (38)- 
(40). This is obtained when, say, an uus state can turn into E + with probability 3/4 
and into A (via E + (1385) — > An + decay) with probability 1/4 0. In the case of Eqs. (41) 
and (42) E~ can not be produced and we assume that E + = 0.3A. 

The values of v , v„ and Vqq 3X6 determined directly by the corresponding coefficients 
of Eqs. (38)- (42), together with the probabilities to fragment a qqs system into E + + E~ 
and into A, given above. For example, we have for incident uu diquark and secondary 
proton : 

Vo = AL 3 ,v q = 3L 2 ,v qq = 2L; (47) 
for incident ud diquark and secondary proton : 

v = AL 3 ,v q = 2L 2 ,v qq = L; (48) 

for incident uu diquark and secondary A : 

12 4 1 

= Yq SIj2 ,Vq = T6 SL ,Vqq = 4 S ' ^ 49 ^ 



for incident ud diquark and secondary A : 

12 qr2 

bL ,v a 

1.6 q 1.6 



12 4 

^0 = -TT SL2 > V q = T^ SL 1 V 1<1 = S ' ( 50 ) 



3 These probabilities are in disagreement with simplest quark statistic rules E8 
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for incident u quark and secondary H : 

v = 3S 2 L ,v q = 0; (51) 
for incident d quark and secondary 2~ : 

v = 3S 2 L , v q = S 2 ; (52) 
and incident SJ and secondary Q~ : 

v = S 3 . (53) 

The probability for a process to have n cutted pomerons was calculated using the 
quasi-eikonal approximation ]1|, 2T| : 



= CTn/ 5]<T», CT n = (1 - e~ Z 51 TiO ' ( 54 ) 

n=l nZ k=0 fc! 



i? 2 + a'f 

with parameters 



: J^ 7 e Ag , cr P = 87r 7 e A? , £ = ln(s/l GeV 2 ) , (55) 



A = 0.139 , a' = 0.21 GeV"' , lpp = 1.77 GeV"' , 7?rp = 1.07 GeV"' , 
R 2 pp = 3.18 GeV" 2 , R 2 np = 2.48 GeV" 2 , C pp = 1.5 , C np = 1.65 . 

The model parameters for quark and diquark distributions and their fragmentation 
are presented in Table 1. They were mainly taken from the description of the data in 
Refs. [HI, HI. 
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Table 1 



The values of the parameters used for the calculations in QGSM. 



Parameter 


value 


a N 


1.8 




0.18 


di 


12 


Mr 


0.5 


a B 


-0.5 
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Figure captions 



Fig. 1. Cylindrical diagram corresponding to the one-pomeron exchange contri- 
bution to elastic pp scattering (a) and its cut which determines the contribution to 
inelastic pp annihilation cross section (b) (string-junction is indicated by a dashed line), 
c) the diagram for elastic pp-scattering with an exchange by SJ in the t-channel and its 
s-channel discontinuity d). 

Fig. 2. Three different possibilities of secondary baryon production in pp interac- 
tions: string junction together with two valence and one sea quark (a), together with 
one valence and two sea quarks (b), together with three sea quarks (c). 

Fig. 3. The spectra of secondary protons (a), and antiprotons (b) in pp collisions 
at 100 and 175 GeV/c [|] and at 400 GeV/c pf (c) and (d) and their description by 
QGSM. Secondary proton (e) and antiproton (f) yields at ISR energies |TT| at 90° in 
c.m.s. and their difference (g) together with QGSM model predictions. In all cases the 
calculations with e = 0.05 are shown by solid curves and the variants with e = 0.2 are 
shown by dashed curves. 

Fig. 4. The spectra of secondary antiprotons in ir~p collisions (a) and of protons in 
n + p collisions (b) at lab. energies 100 and 175 GeV || and its description by QGSM. 
The differences dN/dx oip — p yields in \{ix + p + ix~ p) collisions at 158 GeV/c compared 
to preliminary data from NA49 |K| (c). In all cases the calculations with e = 0.05 are 
shown by solid curves and the variants with e = 0.2 are shown by dashed curves. 

Fig. 5. The spectra of secondary A (a) and A (b) in pp collisions (data are taken 
from 0) and its description by QGSM model. The calculations with e = 0.05 are shown 
by solid curves and a variant with e = 0.2 are shown by dashed curves. 

Fig. 6. The asymmetries of secondary A/A (a), (b), and Q/Q (c), in 7i~p 

collisions at 500 GeV/c |L2| and its description by QGSM model. In all cases the calcu- 
lations with e = 0.05 are shown by solid curves and the variants with e = 0.2 are shown 
by dashed curves. 
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